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Abstract 

We develop a first-principles model of thermally-activated cross-slip in magne- 
sium in the presence of a random solute distribution. Electronic structure meth- 
ods provide data for the interaction of solutes with prismatic dislocation cores and 
basal dislocation cores. Direct calculations of interaction energies are possible 
for solutes — K, Na, and Sc — that lower the Mg prismatic stacking fault energy to 
improve formability. To connect to thermally activated cross-slip, we build a sta- 
tistical model for the distribution of activation energies for double kink nucleation, 
barriers for kink migration, and roughness of the energy landscape to be overcome 
by an athermal stress. These distributions are calculated numerically for a range 
of concentrations, as well as alternate approximate analytic expressions for the 
dilute limit. The analytic distributions provide a simplified model for the max- 
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imum cross-slip softening for a solute as a function of temperature. The direct 
interaction calculations predict lowered forming temperatures for Mg-0.7at.%Sc, 
Mg-0.4at.%K, and Mg-0.6at.%Na of approximately 250°C. 
Keywords: magnesium alloys; dislocations; plastic deformation; cross-slip; 
density functional theory 

1. Introduction 

Increased interest in the light-weight structural metal magnesium[[T]| to re- 
place aluminum or steels in automotive applications^ has focused attention on 
a variety of metallurgical issues, including formability. Current Mg alloys re- 
quire temperatures near 300°C for forming to activate the five independent slip 
systems required by the von Mises criterionLSJ; this is in part due to the large 
anisotropy between basal and prismatic slipflH. Cross-slip of a-type dislocations 
from the easy (0001) basal plane onto the hard (0110) prismatic plane requires 
large stresses or high temperatures. Experimentally, few solutes have been found 
to lower the stress for cross-slip: Al and Zn lower the stress at low (below room) 
temperatures [5], while Li can lower the cross-slip stress in both regimes^ |71 
[HI |9l. The difficulty of performing experiments to measure cross-slip stresses 
for alloys — ^requiring single-crystal samples oriented for prismatic slip — is com- 
pounded by the possibility that, like solid-solution softening in BCC alloys [fTOll . 
it may occur over a limited concentration and temperature range. Hence, new 
state-of-the-art first-principles prediction of solute/dislocation interactions cou- 
pled with predictive computational modeling of thermally-activated cross-slip in 
the presence of solutes is necessary to guide the design of new alloys. 

Couret and Caillard in situ experimental measurements[[IIl [111 found that at 
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high (above room) temperatures, cross-slip in magnesium is the result of a double- 
kink nucleation (also called "jog-pair"[fT3l) mechanism. A basal screw dislocation 
constricts at kinks of height c and spreads on two neighboring basal planes [fT4l . 
At low (below room) temperatures, cross-slip instead occurs by constriction and 
bowing of a screw dislocation — the Friedel-Escaig mechanism [[T5l [T6ll : as we are 
interested in the problem of forming near 300°C, we consider only the double- 
kink nucleation mechanism. We model double-kink nucleation by computing 
the geometry and formation energy of a single basal-to-prismatic kink with a 
validated embedded-atom potential [fTTl : together with the Peierls stress, we ac- 
curately reproduce the experimental cross-slip stress over the 300K-700K tem- 
perature range. The starting geometry and energetics are combined with first- 
principles modeling of the changes in energy from basal to prismatic cores due 
to substitutional solutes in magnesium. Finally, we compute — numerically and 
with an analytic approximation — the distribution of double-kink activation bar- 
riers and energy barriers for kink migration to predict the stress for thermally- 
activated cross-slip with solute concentration and temperature for K, Na, and Sc. 
In order to do direct substitution of solutes, we can only consider solutes which 
do not increase the energy of the prismatic core. Finally, we use our model with 
solute data for softening to predict the maximum possible reduction in forming 
temperature for a binary alloy. 

2. Computational methodology 

To accurately compute the interaction of solutes with Mg dislocation cores, 
we use density-functional theory with flexible boundary conditions for a stress- 
free dislocation core. Density functional theory calculations are performed with 
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VASp ffT^fT^ . using a plane- wave basis with the projector augmented- wave (PAW) 
method [|20| and potentials generated by Kresse[[2TI|. The many-electron exchange 
and correlation effects are described by a generalized gradient approximation [|22]| . 
and a plane- wave kinetic energy cutoff of 273eV ensures accurate treatment of the 
potentials. The electronic configurations of the PAW potentials are Mg: [Ne]35^, 
K: [Mg]3p''As\ Na: [Ne]35\ Sc: [Ar]4s^3d^. The dislocation calculation use 
a 1 X 1 X 16 /:-point mesh with a Methfessel-Paxton smearing of 0.5eV to give 
an energy accuracy of 5meV for bulk Mg. The dislocation geometry and flex- 
ible boundary condition method ll23l l24l l25l use a similar geometry to previous 
calculations for basal dislocations [|26ll with a reshaped region I for prismatic split- 
ting; we relax atomic positions until the forces are less than 5meV/A, where the 
lattice Green function is used to displace atoms in the outer regions. To pro- 
duce a prismatic screw dislocation, we begin relaxation from the anisotropic dis- 
placement field from two b/2 partial screw dislocations (for total Burgers vector 
b = |[21 10]) separated by c in the prismatic plane. 

Determining the geometry and formation energy for a basal to prismatic kink 
requires a validated classical potential [fTTl. This EAM potential, optimized for liq- 
uid/solid interfaces, also accurately reproduces stacking fault energies and dislo- 
cation core geometries compared with density-functional theory calculations [|27l 
l28ll . The geometry of a Z? = |[2110] kink of height c is computed from a pe- 
riodic cylinder with axis A^|[2110] -I- c[0001]; this produces a periodic array of 
kinks with density A^"^. The inner cylinder radius is 22Z7 plus an outer layer of 
5b fixed to the initial positions from anisotropic elasticity for a mixed dislocation. 
The initial geometry relaxes to a metastable mixed dislocation configuration. To 
produce a kink, we (1) apply a small Escaig shear strain of -0.006 between the 
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I [21 10] and a[0110] directions to promote basal spreading, (2) relax, (3) reverse 
the strain, and (4) relax. The energy of the configurations varies as a linear func- 
tion of N: constant kink energy plus the line energy with N. The energies of the 
core cylinder at = 60, 80, 100, 120 go as (0.515 + 0.9 lOA^) eV, and we identify 
the single kink formation energy Ef = 0.515eV. The core cylinder — where the to- 
tal energy becomes linear in radius — is I2.5b; if instead the sum is carried out to 
22b, the formation energy is 0.510eV. The EAM potential has a zero-temperature 
prismatic Peierls stress t* = 140MPa; this is the minimum stress required for a 
straight a-type screw dislocation to glide in the prismatic plane. We validate with 



the experimental data for cross-slip from 300-7 OOK in Section 3.3 



3. Results 

3.1. Prismatic screw dislocation core 

Figure |3] shows the core of a prismatic screw dislocation from first-principles. 
In Mg, the a-type screw dislocation splits into partials in the (0001) plane; split- 
ting onto the (0110) plane is a metastable higher energy configuration. The basal 
stacking fault energy is 34 mJ/m^ compared with the prismatic stacking fault 
energy of 218 mJ/m^; the metastable point for the prismatic stacking fault is 
|[1120] -I- 0.065c[0001]. We identify 5c = 8 sites in the core with large screw 
dislocation content where solutes may substitute and lower (or raise) the disloca- 
tion core energy. Unlike the basal screw dislocation, there is very little volumetric 
strain and the dislocation core spreads rather than splitting into partials. Because 
the prismatic core is metastable, solutes can destabilize the core; in that case, re- 
laxation of the geometry will split the dislocation core on the basal plane. Here 
we use direct calculation of interaction energy for attractive solutes. 
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3.2. Direct solute/dislocation interaction 

Figure |4] shows the difference in energy for solutes substituted at sites in a pris- 
matic dislocation core relative to substitution in a basal dislocation core. Three 
solutes decrease the prismatic stacking fault energy in Mg and produce attrac- 
tive dislocation interactions: K, Na, and Sc. Note that repulsive interactions are 
still possible as the sane site in the basal dislocation core can be more attractive 
than in the prismatic dislocation core. Of these three solutes, K has the strongest 
interaction in both dislocation cores due to the large size misfit and change in 
stacking fault energies ("chemical" misfit [|26il ). The statistics of the different site 
energies — mean site energy over the eight core sites AE and standard deviation 
6E — are included in Table [T] These parameters enter into the analytic approxima- 
tion to the statistical model we derive, and from which we compute the dimen- 



sionless softening parameter Pdk in Section 3.4, Tath/T* in Section 3.5 and their 
ratio X- 

3.3. Kink geometry and enthalpy 

Figure [T] shows the relaxed kink geometry. As kinks are "decorated" with 
solutes, we need to know the portion of the kink that is prismatic; these sites 
will have their energies changed by the energy difference between a prismatic 
and basal core containing a solute. The kink has a height of c, and takes 30b to 
step from one basal plane to a neighboring basal plane; over that length, 15b is 
a prismatic core; hence, for ^^ink = 15, we have S = Sc^idnk = 120 sites which 
can be occupied by solute atoms. We consider kinks with random distributions 
of solute, and derive the probability of having a kink with energy E; we then 
determine the distribution of double -kink nucleation energies and kink-migration 
barriers as well as the athermal stress required for kink mobility. 



Figure |2] compares the experimentally determined cross-slip stress for pure 



Mg ll29]| with the calculated double-kink nucleation enthalpy. In a manner similar 
to Kocks et g/. OOll . we expect the double-kink nucleation enthalpy to follow 



where Ef is the formation energy of a single kink, and r* is the prismatic Peierls 
stress; the exponent of 1/2 corresponds to the elastic interaction of two well- 
separated kinks [|3T]| . To validate this expression, we use the Orowan equation 
to relate the enthalpy to the plastic strain rate s for prismatic shear 



where is the dislocation density (taken as lO^cm"^), and Vdk is the double kink 
nucleation attempt frequency. The double -kink nucleation (or jog-pair 111 311 ) model 
for cross-slip is valid in the "high" (above room temperature) range IfTTl [T2l [T4ll 
where the stress is low enough that constriction on the prismatic plane is limited 
to the kinks. At low temperatures, cross-slip occurs through bowing and this 
treatment does not apply. The attempt frequency is difficult to compute accurately; 
however, it can be estimated as a typical phonon frequency (~ lOTHz) divided 
by ~ 45 (the number of atoms in two kinks), or instead fit to the experimental 
data in the 300-700K range. A single parameter fit with the experimental strain 
rate e = 1.66 x lO'^^s"^ gives Vdk = 15.4GHz, which is remarkably close to our 
simple estimate counting degrees of freedom. With this one parameter, we are 
able to reproduce the experimental cross-slip stress from room temperature up to 
~ 400°C, 




(1) 




(2) 



T(r) = T* 1- 




(3) 
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where 




2£f 



= 62.5meV 



(4) 



In(capxVdk) - Ins 



for our dislocation density and strain rate. For notational convenience, we intro- 
duce the parameter ordk as 



for pure Mg, a^^^ = k^T/H'^^ from Eqn. (jsj). The excellent agreement with both 
tensile tests and in situ experiments validates our use of a double-kink nucle- 
ation model for cross-slip above room temperature; it should be noted that there is 
marked deviation going to absolute zero, as the mechanism for cross-slip changes. 

3.4. Distribution of double-kink nucleation energies 

Each nucleated pair of kinks in a solute field requires a total energy equal to 
the formation energies plus the change due to the presence of solutes in the kink; 
the distribution can be computed numerically. Each kink has 5 = 120 sites in 
the prismatic core of the kink that may be occupied by solute atoms; we assume 
that the total energy change for the kink is the sum of all the individual energy 
changes, which are given by the site occupancy (either or 1) multiplied by the 
energy of that site. We assume translational invariance down the length of the 
kink, so that there are only 5 c = 8 unique site energies to consider labeled AE',- 
for z = 1 . . .5c — these are the site energy changes in Figure |4j Then, in the kink 
core, each row of sites has solutes (between and i^^k = 15), which contribute 
energy niAEt; the total occupancy of the S core sites is n = 2; and the energy 
is £■ = Yiif^i^Ei. Hence, we can write the number of possible configurations 
involving no out of S kink core sites occupied with energy change Eo as 




(5) 




(6) 
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where 6(Eo - E) is the Dirac delta function, 5„ „q is the Kronecker delta, and the 
final term accounts for the multiplicity of occupancies along each row. For numer- 
ical convenience, we approximate the delta function with a smoothed Gaussian 
with width lOmeV. From Eqn. (|6]), the fraction of double kinks that can form with 
energy lEf + £ for a random solute distribution concentration Cs is 

G''\2Ef + E,c,) = J]J] ^""(1 - c,)^^-"-'" dE'giE', n)g(E - E', m). (7) 

«=0 m=0 

Figure |5] shows the numerical double-kink nucleation energy distributions for sev- 
eral concentrations. The average nucleation rate for double kinks at stress r and 
temperature T is 

dEG^\lEi + E,c,)QX^ 

CO 

assuming C^^ilEf + E, Cg) « for £ < -lEf, otherwise the integral must be split 
at £ = -2Ef. 

An analytic approximation can be derived for the distribution of energies and 
the average nucleation time. Eqn. ([6]) (and hence, Eqn. (|7])) can be alternately 
viewed as the distribution of the sum of n random variables, where each variable 
is the site energy. From the central limit theorem, this distribution will be normal 
with mean value of nAE and standard deviation y/nSE (for average interaction AE 
and standard deviation 6E;c.f. Table [T]), 

Ct(2E,..,.) = |:(J<(l-.,-— e.p 

This form matches the appearance of distributions in Figure [5] The analytic ex- 
pression of Eqn. ^ can integrated in Eqn. ([s]), using = (k^T)'^ and from 



2E{ + E 



(8) 



(E - nAEf 
2n6E^ 



(9) 
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Eqn. Q 



Vdk 



exp 



dE- 



-fc(2£f + E)- 



(E~nAEf 
InSE^ 



n=0 



(10) 



9 c 



= Vdke {l 

where the first term is the nucleation rate in the absence of solutes, and the term 
in braces is the change in rate due to solutes. 

The analytic rate equation can be simplified to give the softening stress with 
concentration in the low concentration limit. The change in cross-slip stress can 
be understood as a change in the stress necessary to have the same nucleation rate 
at a given temperature in the absence of solutes; that is, the cross-slip stress r has 
a corresponding a^^ such that 



C +C g-/'"dkA£g08adk'5£)'/2 



where = k^T jH'^^^. For small Cg, dadk/dc^ is 



'dk 



^{an,-a°^^2Ef/{2S) _ ^ _ ^^^~/ian,AE^(fiaii,SEf/2 

1 + ^(adk - O - 1 - c. + c,.-^'^Sk^.^<^^)^/^ 



Mk^ 



dc, I32E 



exp 



' AE ^1 I6E\ 
< ^1 2 \H°J , 



2\ 



Finally, as r = r*(l - Cdk)^, the change in cross-slip stress is 



— = -T 2(1 -Q-dk)-^ 
acs dCs 



2Ef 



45 



exp 



' AE ^1 I6E\ 



2\ 



dk^"^dk ■ -f^dk 



(11) 



(12) 



(13) 
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where P^k is the unitless "softening potency" (c.f. Table [T]). Note that P^t > 



soften cross-slip provided sufficient attractive sites are available. 

3.5. Distribution of kink migration energy barriers 

Solutes "roughen" the energy landscape for kinks and provide two barriers 
to the motion of kinks: a minimum (athermal) stress required for kinks to mi- 
grate preferentially down the dislocation line, and the energy barrier over the 
length of a kink[32]. Solutes provide local changes in energy as a kink moves 
over a single Burgers vector; a minimum stress — r^(Cs) — is necessary to over- 
come this short-range change in energy. This is given by the average roughness 
of the energy landscape: the energy changes as 5 c sites with solutes "leave" the 
kink and another 5 c sites "enter" the kink; the energy change for solutes entering 
is +AEj, and is -AE, for those leaving. If the width of the solute interaction is 
approximated as a Gaussian with area 4b^, the stress to overcome should go as 
^/2/e ■ (energy difference)/4Z?^. Analytically, we consider a distribution of energy 
differences as the sum of n random variables for the energies; however, both +AEi 
are equally likely so the mean is and standard deviation 16E + AE 1 . The av- 
erage energy difference is symmetric around n = Sc', hence, the standard deviation 

/ 2x1/2 

of our normal distribution for n = 0..5c is y/n{6E + AE I and for n = 5c. .25c it 



is V25 c - n(6E^ + AE ) . From this, we can approximate the average absolute 



requires AE < 



(6E/H^^^)^, so that solutes with repulsive interactions can 
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\AE\ = 



energy change \AE\ as the kink moves by one lattice spacing, 

fl I o 2x1/2 

J -(IS ,)c,[6E^ + AE ) 



(14) 



and then the athermal stress is 



2\AE\ 2 25, 



/ 2x1/2 

[6E^ + AE) c, = t.. 



ithCs 



(15) 



e 4P 4b^ 

where the last expression is a simple analytic model in the limit of low concen- 
tration. The numerical distribution of energies is the athermal distribution in Fig- 
ure [5] In addition to this athermal stress, there are energy barriers as a kink moves 
by its length ^tink due to solute occupancy changes. For a random distribution of 
solutes at concentration Cs, the fraction of energy barriers E is 

S S ^oo 

G^'^iE, Cs) = J]Yj ^"^"^^ ~ I dE'giE', n)g{E + E\ m). (16) 

«=0 m=0 

Note the sign change for g{E + E',m) compared with Eqn. (|7]). The enthalpy 
barrier to escape an energy well £ > is (assuming a Gaussian of width 2£]ankb 
and kink height h = c) 



1 - 



2i^^Jy'hT-T\c,) 



-ill 



(17) 



for T > T^(Cs), from [|32l. Define 



E, = (T-T\c,)y 



2£khikb^ 



We 



(18) 
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then the time needed to overcome all barriers along the dislocation line is 



dEG^'^{E,c,)+ I J£G™(£, Cs) exp 



JEr 



km/- 



1-^ 

E 



(19) 



Note that is finite only when r > t^(Cs). 

There is a minimum solute concentration — and hence, athermal stress — necessary 



for thermally- activated kink migration to affect the cross-slip stress. Eqn. (16) is 
a distribution of the sum of n random energies; however, as solutes enter and 
leave the kink, our possible energies are iAE, for a distribution with mean and 
Standard deviation \ 6E + AE 1 as with the athermal barrier. Due to symmetry 
around n = S , the standard deviation of our normal distribution for no = 0..S is 
^/n[6E^ + AE^'^ and for n = 5. .25 it is V25 - n{6E^ + AE^'^ . From this, we 
approximate G^{E, Cs) as normal with mean and variance 



(E )km - 



«=0 ^ " ^ 71=5 + 1 ^ " ^ J 

x(2S)c,(5E^ + AE^) 

(20) 

where the approximation has less than 1% error for Cs < 45% with S = 120. 

1 / 9 2\l/2 1 /-^ 

Hence, the standard deviation is ^ ylSiSE + AE 1 . With this approxi- 
mation, the stress necessary for thermally activated kink migration to require a 
larger stress than athermal kink migration is when at r = t^(Cs) is slower than 
required by the Orowan equation: s > capj^t^^ (c.f. Eqn. p])). At the athermal 
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stress, there is no reduction in enthalpy for E > 0, and so 

1 r°° 

?ionVdk=2+ dEG^'^{E,c,)^'' 



1 r 



dE ^- ^ ^ (21) 



exp 



U\2S)(5E^ + ^Ey, 



r 

where the first approximation is the use of a normal distribution and the second is 
valid when the exponential term is larger than 1 . Then, the minimum concentra- 
tion c^'" is such that s = capj^t^^; above this concentration, thermally activated 
kink migration will be required. Then 

\/S\2S){5E' + AE'yr = In (^^) 

(hTf 2Ef (22) 



mm 

*^ c 



(6E^ + AE^) ^dk^ 

We rewrite this in terms of the minimum value of athermal stress, as t^(Cs) « 

^.^_(kBTf(6E^ + AE) 2E, 

''"min * "^athC" " / , 3 ^ 775" " ^^^^ 

If T^(Cs) < rj^jj^, then thermally activated kink-migration will lead to further strength- 
ening. For the solutes we consider here (c.f. Table [T| and Figure [6]), t^j^^ is > 
15MPa at 300K and > 40MPa at 600K; hence, thermally-activated kink- migration 
only limits softening outside of the stress and temperature range of interest. Above 



c™", thermally-activated kink-migration controls cross-slip. The integral in Eqn. ( 19 1 
can be computed in closed form with the approximation 



G'^(£,c,) « exp(-£V(2<£')km))/V2;r<£2) 



km 
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from Eqn. ( [20] ) and with 

(2 - V3) (^) (3 + V3)£J2 



E 



(24) 

:£>(3+V3)£J2 



The closed-form expression for Eqn. ( 19 1 is unwieldy, and omitted for clarity. 



3.6. Prediction of cross-slip stress with concentration and temperature 

Figure |6] shows the predicted cross-slip stress with concentration and tem- 
perature predicted by our numerical solution of the Orowan equation, and ana- 
lytic approximations for softening and athermal hardening. The Orowan equation, 
Eqn. (|2]), is modified for the two thermally-activated cross-slip processes as 



e = 



(^)z.vxfek+ur- (25) 



This relates the plastic strain rate to the time to nucleate a pair of double-kinks 
and migrate the length of the dislocation line. The analytic curves are for double- 
kink nucleation and athermal hardening, both of which are approximately linear. 
For all solutes, the analytic approximations are reasonable at lower Cs, with more 
significant deviations at higher concentrations. All three solutes show softening 
for low concentrations which becomes hardening at higher concentrations; the 
range of solute concentration that leads to softening decreases at higher temper- 
atures. These predictions suggest possible alloying concentrations that can lead 
to lower stress for thermally-activated cross-slip, and hence decrease the plastic 
anisotropy. At 600K, the cross-slip stress is 4.1MPa from our model; that cross- 
slip stress occurs at 545K for Mg-0.4at.%K, 560K for Mg-0.6at.%Na, and 565K 
for Mg-0.7at.%Sc. 
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We derive approximate analytic expressions for the dilute concentration limit 
for alloy design. The optimal solute concentration occurs when the minimum 
stress for double-kink nucleation matches the athermal stress; hence, 



T*(l - <)(1 - < - <PdkC°) = TathC° 

(l-<f-(l-<X/',,c« = I^'^-° 

T 



(26) 



(1 - <f 



(1 -<)<^dk + Tath/T* 



and then the cross-slip stress is TathCj . Finally, the cross-slip stress for pure Mg at 
one temperature can now be achieved at a minimum "equivalent temperature" in 
an alloy. The pure Mg forming temperature T^orm is written as a^,^ = ^BT^form/^dk' 
and the minimum equivalent temperature rminimum as adk = ^BT'minimum/^di,; then. 



T*(l-0 



) — 

(1 - Qrdk)ffdk^dk + Tath/r* 



1 - adk 

1 - a" 
^ "dk 



2 



- 1 



— ^^(1 - Q'dk)ffdk- 

Tath/T* 



(27) 



Define = /'dk/(Tath/T*), the unitless ratio of softening potency to athermal hard- 
ening; the higher this parameter, the more the forming temperature of the alloy 
can be lowered. The solution to the quadratic formula is 

l-adk= ^-^ (28) 

1+^(1 -<)2 

This is plotted in Figure |7] for a variety of;^^ values. Comparing to our numerical 
values, the analytic approximation gives equivalent 600K forming temperatures 
of 530K for Mg-K, 543K for Mg-Na, and 560K for Mg-Sc. Note that for Tform = 
600K and rminimum = 300K (room temperature forming) requires^ = 43.3, which 
would be an attractive interaction to a kink of at least ~ 4//jj^ = 250meV. 
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4. Conclusions 

Our first-principles calculations of solute interaction with prismatic and basal 
dislocations combined with our computational model of basal to prismatic cross- 
slip above room temperature allows us to make the first prediction of cross-slip 
softening of Mg in three binary alloys with lower forming temperatures. This 
model accounts for the random distribution of solutes in a dislocation core with 
different possible energies, and the direct calculations can be used above the di- 
lute concentration limit. At the same time, an approximate analytic model con- 
nects average solute interaction energies to the dilute limit for near-quantitative 
predictions of cross-slip softening, and to determine the required interaction en- 
ergies that lead to softening. This predicts possible formable Mg alloys and sets 
the stage for future work covering a wider range of solutes; moreover, the aspects 
of the approach can be applied to other thermally-activated plastic deformation 
mechanisms. Replacing direct calculations of interaction energies with a geo- 
metric model connecting the prismatic dislocation core to changes in prismatic 
stacking fault energy could also treat other substitutional solutes. 
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Figure 1: Kink geometry from Sun potential. The atom colors indicates changes in atomic energy 
relative to bulk Mg. 

Table 1: Direct solute/kink interaction energy statistics and parameters for analytic approxima- 
tions. The analytic softening and hardening depend on the distribution of solute energies around a 
kink; from these, the (unitless) solute softening potency Pdk and (unitless) athermal slip prefactor 
Tath/"!"* are derived. The ratio of the two factors, j^', determines the maximal amount of softening 
that is possible in a given alloy; larger values indicates more potential softening (c.f. Figure |7]l. 
The linear analytic approximations are reasonable for Cj < 2% (c.f. Figure [6]l. 

interaction [meV] strength param. 

K -45.8 44.0 63.5 48.5 6.67 7.27 

Na -37.0 22.3 43.2 26.9 4.54 5.93 

Sc -19.6 38.1 42.8 18.9 4.50 4.20 
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Figure 2: Cross-slip stress with temperature from experiment and double-kink nucleation enthalpy 
value of Eqn. (|3]l. The single kink formation energy and prismatic Peierls stress are computed 
directly from the Mg EAM potential; the Orowan equation determines a cross-slip stress for a 
given dislocation density and plastic strain rate. There is a single unknown parameter Vdk which is 
determined by a least-squares fit to the experimental data; simple counting of degrees of freedom 
suggests an order of magnitude of ~ 10 '"Hz as found in our fit. 
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Figure 3: Basal and prismatic a-type screw dislocation core geometries. The differential displace- 
ment plots (arrows) indicate relative displacement of rows by up to bjA to show the small amount 
of splitting into partials. The coloring is a linear-interpolation of Nye tensor density indicating up 
to a total displacement of ft/2. The metastable prismatic core has a small amount of spreading in 
the prismatic plane. 
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Figure 4: The relative energy of a K, Na, and Sc solute in a prismatic dislocation core compared 
with the same site in a basal dislocation core. The change in energy can promote (negative en- 
ergies) or prevent (positive energies) the formation of a double kink. The interaction energy of a 
single solute with a single dislocation enters into the distribution of double-kink formation ener- 
gies, kink-migration barrier energies, and athermal kink migration stress. 
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Figure 5: Probability distributions for all three processes (double-kink nucleation, thermally- 
activated kink migration, and athermal kink migration) and solutes (K, Na, Sc). The probability 
distributions are determined numerically from the solute interaction energies assuming random 
solute distributions. For "large" solute Cs (compared to the inverse of the number of available kink 
size 5 '), the distributions become close to normal, and approximate analytic expressions can be 
derived. 
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Figure 6: Cross-slip stress with concentration for K, Na, and Sc between 300-700K. Increasing 
temperature lowers the pure Mg cross-slip stress; small additions of attractive solutes increase the 
double-kink nucleation rate, leading to softening. The lower bound for all curves is the athermal 
stress T^(cs); only at higher concentrations does thermally-activated kink-migration need to be 
activated to achieve the required plastic strain rate. The dashed lines show the linear softening 
analytic approximation, and the dashed-dotted curves are the linear analytic approximation to the 
athermal stress. All cases show a minimum cross-slip stress for a given temperature, and the 
concentration range for cross-slip softening decreases with increasing concentration. This allows 
us to identify an "equivalent temperature": the lowest temperature where the binary alloy has the 
same cross-slip stress as pure Mg has at a higher temperature. Hence, forming of these binary 
alloys can be done at a lower temperature, provided the solute concentration is properly tuned. 
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Figure 7: Analytic equivalent minimum temperature for different solute interaction parameters 
X - Ptik/iTmh/T*). Contours for each;^' are shown, increasing logarithmically. In the analytic ap- 
proximation, the ratio of double-kink softening to athermal hardening determines how effective a 
solute can be at reducing the cross-slip stress for a Mg alloy. The forming temperature for pure Mg 
defines the necessary cross-slip stress; the minimum temperature that a solute can give the same 
cross-slip stress is plotted against that forming temperature for a range of parameters. The inset 
shows the decrease in forming temperature relative to 600K. Note: the thermally-activated cross- 
sUp model is only valid for temperatures above room temperature, and the analytic approximations 
are accurate for Cs < 2%. 
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